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Abstract 



We consider a two-parameter (c,c) family of quantum integrable Hamiltonians for a 
chain of alternating spins of spin s = 1/2 and s = 1. We determine the thermodynamics 
for low-temperature T and small external magnetic field H, with T « H. In the antifer- 
romagnetic (c > , c > 0) case, the model has two gapless excitations. In particular, for 
c = c, the model is conformally invariant and has central charge c v i r = 2. When one of 
these parameters is zero, the Bethe Ansatz equations admit an infinite number of solutions 
with lowest energy. 
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The one-dimensional Heisenberg model, like the hydrogen atom, has served as the 
guiding example for a very large body of both experimental and theoretical work. Progress 
has recently been made on closely related models, consisting of chains with alternating 
spins, such as spin 1/2 and spin 1. On the experimental side, materials (e.g. 1 , [MnCp^] 
[TCNE]) have been synthesized which, at temperatures above a certain transition temper- 
ature T c , behave as one-dimensional ferromagnets of alternating spins. On the theoreti- 
cal side, quantum integrable models of chains with alternating spins have recently been 
constructed 2 . In this Letter we investigate the thermodynamics of a two-parameter fam- 
ily of such integrable models. Depending on the values of the parameters, we find both 
antiferromagnetic and ferromagnetic behavior. When one of these parameters is zero, the 
Bethe Ansatz equations admit an infinite number of solutions with lowest energy. 

We consider a system of N spins |c?2 , |<74 , • • • , |oW of spin 1/2 and iV spins s i , s 3 , 
• • • , S2N-1 of spin 1 in an external magnetic field H(> 0) with the Hamiltonian Ti given 
by 

H = cH + cH- HS Z , (1) 

where S z = E^i + £n=i s 2n-n 
1 N 

H = - g ( 2 ° 2n ' * 2n+1 + ^ ( 2 ^n+2 • S 2n +1 + 3) , (2) 

n=l 

1 N 

jj = " ^ (2<?2n • hn-l + 1) [(1 + S2n-1 ' «2n+l) (2a 2 n • «2n+l + 1) + 2] , (3) 

y — * 

n=l 

and c and c are real constant parameters. (In this paper, bars and tildes are interchanged 
with respect to Ref. 2.) Note that the Hamiltonian contains both nearest and next-to- 
nearest neighbor interactions. We assume periodic boundary conditions: om = &2n+2N 

and S 2n +1 = S2n+1+2AT- 

The corresponding energy eigenvalues are given by 2 

E = cE + cE-H QiV - M^j , (4) 

2 



where 

where the variables Xj satisfy the Bethe Ansatz (BA) equations 



N M 

^±j^±i = Aj-Afc + i J = 1 ,..., M . (6) 

We consider here a strictly alternating arrangement of spins, with spins 1/2 at even sites 
and spins 1 at odd sites. For any other ordering of the spins, one can construct a corre- 
sponding Hamiltonian which has the same energy eigenvalues and BA equations. 

This system of equations admits the same string solutions that are found for the 
Heisenberg model. In the thermodynamic limit, the model is characterized by particle 
densities p n (X) and hole densities p n (A). Following the standard procedure (see, e.g., Refs. 
3, 4), we find that these densities obey the constraints 

cso min(n,2) 

a n + a "+3-2Z; ( 7 ) 

m=l 1=1 

where 

min(n,m) — l 

Anm(X) = 5 nm 5(X) + (1 - 5 nrn ) a| n _ m |(A) + a n+m (A) + 2 ^ a| n _ m | +2 z(A) , (9) 

i=i 

and * denotes the convolution (/ * g) (A) = dX' f(X — X')g(X') . 

The thermodynamic Bethe Ansatz (TBA) equations describing the equilibrium system 
at temperature T are 

oo / min(n,2) \ 

Tln(l + e e "/ T ) = ]T [A nm * T In (\ + e~^' T ^ - 2tt ( ca n + c ^ a n+3 - 2l \+nH, 

(10) 
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m=l \ 1=1 



where 

e n (A) = Tln(p n (A)/p n (A)) . (11) 

(The particle and hole densities are understood to be those at equilibrium.) The equilib- 
rium free energy is given by 

F = 2Ne -NT f dX s(A) In |Yl + e £l(A)/T ) (l + e £2(A)/T )j , (12) 

J — oo 

where s(A) = l/(2ch7rA) , and 2e = c [| - (21n2 + 4-7r)] +c [-\ - (6 - tt)] . 

In the high-temperature limit with zero field, these equations give the expected value 
for the entropy, namely 5 = iVln6. We now consider small values of T and H, with 
T « H. We define e n (X) = lim^^o e n (A) . Keeping in mind that e n > for H > and 
n > 3, we obtain from the TBA equations the following system of linear integral equations 
for Si and e^'- 

E\ = —27TCS + s * e 2 , 

e 2 = -H-2ncs + h * e£ + s * e+ , (13) 

where h = s * a\ , and the superscript + on a function denotes the positive part of that 
function; i.e., e + = \ (e + \e\) . 

The qualitative behavior of the solutions depends on the sign of the parameters c 
and c, and the various cases must be studied individually. Let us first consider the anti- 
ferromagnetic case c > 0, c > 0. In this case, for H = the solutions £i(A) and 62(A) 
are readily found, and are seen to be negative for all A. Hence, for the ground state, 
Pi (A) = /02(A) = s(A) , and all other particle and hole densities are equal to zero 5 . Thus, 
the ground state is a "sea" of strings of length 1 and a "sea" of strings of length 2, in agree- 
ment with the alternative analysis of Ref. 2. This corresponds to the antiferromagnetic 
ground states of integrable chains 4 of spin s = 1/2 and s = 1. 
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For H small but nonzero, we can solve for e\ and e 2 by generating from (13) a system of 
Wiener-Hopf equations. To this end, we define a n to be the zeros of £ n (A), i.e., e n (a n ) = , 
for n = 1 , 2 . We assume 



1 

TV 



In H + In K n + O 



1 



In if 



(14) 



where the constants K n (which are independent of H) have still to be determined. Intro- 
ducing the functions 

e™* K n e n (X + a n ) A>0 ^ 



5„(A) t Q A < ' 

we obtain (for H — > 0) the following system of Wiener-Hopf equations 

/•OO 

Si (A) = -2 7 rcKie- 7rA + / dA' s(A - A' + a x - a 2 )S 2 (V) , 

S 2 (A) = - - 2 7 rc« 2 e- 7rA + / dA' [s(A - A' + a 2 - ai)Si(A') + /i(A - A / )S 2 (A / )] , 
2 Jo 

A>0. (16) 

In order to find the leading-order temperature dependence of the free energy, we must 
compute the leading correction to the solutions e n = e n of the linearized equations (13). 
Hence, we set e n (A) = e n (A) + rj n (\) in the TBA equations and expand to leading order 
in r] n , as is explained in Refs. 6 and 7. In terms of the functions 



T n (X) = \ ^-Vn(X + a n ) A>0 
v 1 \ A < ' 

we obtain (for H — > 0) a second system of Wiener-Hopf equations: 

s(A + ai - a 2 ) 



(17) 



S 2 (0) 



poo 

+ / dX' s{X-X' + ai -a 2 )T 2 {X') 
Jo 



r 2 (A) 



/i(A) s(A + a 2 — «i) 
Sp) + ^(0) 



+ / dX' [s(X -X' + a 2 - ai)Ti(A') + /i(A - A , )T 2 (A / )] , 
Jo 



A > 0, 



(18) 



where S' n (0) = dS n /dX 



A=0+ 



Both systems (16) and (18) involve the same 2x2 matrix kernel, which in Fourier 
space is given by 

/ p -iui(a 1 -a 2 ) \ 

^)=( e M-,-l) SM ' US) 

(Our convention for Fourier transforms is that f(u) = f_ 00 d\ e tXuj f(X).) Results from 
Ref. 8 imply that the following factorization exists 

(l-K(u>f) 1 = G+(u) G-(u) , -oo<w<oo, (20) 

where G + (u) and G+ (uj) are analytic in the upper-half complex uj plane with G + (+oo) = 
1, and (for uj in the lower- half plane) G-(uj) = G + (—uj) t . 

Using standard Wiener- Hopf methods, we conclude that the solutions of Eqs. (16) 
and (18) are given (in matrix notation) by 

S{uj) = \ (-L- - -LA G + (,)G-(0) (?) , (21) 



and 



where 



T( U ) = (G+( U ) - 1) [ Ifsm) ' (22) 



S'(0) = - lim u,*S{u) = ^G-{0)( ) . (23) 

|w|— »oo Z V / 

Moreover, the parameters K n introduced in Eq. (14) are given by 

-ftHOo* ^) g -(-)- ig -(°)(;)' (24) 

In order to calculate the free energy per site f = F/ 2iV, we substitute e n = e n 4- rj n 
into the expression (12) for the free energy, as is explained in Refs. 6 and 7. We obtain 

7T 2 T 2 

f = e -H 2 A-—B, (25) 

with 

_ . _ /i / Q'fnWI 

(26) 



A = k t S(m) , B = k t 
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For c 7^ c, the quantities A and B cannot be explicitly evaluated without having 
explicit expressions for the factors G + (u) and G-(u), which we have not yet found. 

For the special case c = c = c, the quantities A and B can be readily evaluated, and 
we conclude (for T « H ) 

It follows that the magnetic susceptibility and specific heat, to lowest order, are given by 



d 2 f 



X = 



1 „ 



C H = —T- 



2n 2 c ' dT 2 



= ^~T, (28) 
h 3c 



OH 2 
respectively. 

This model has 2 ' 9 two gapless excitations, with corresponding speeds of sound v = 2nc 
and v — 2itc. Evidently, the case c = c is the unique case for which the two speeds of sound 
coincide, and the model is conformally invariant. For a critical chain, the low-temperature 
free energy per site is given by 10 ' 11 

/ = e -^T 2 + ---, (29) 

where c v i r is the central charge of the Virasoro algebra and v s is the speed of sound. 
Therefore, from (27) we see that c„; r = 2. Presumably it is no coincidence that precisely 
for the conformally invariant case, explicit expressions for G+(u) and G-(u) are not needed 
to evaluate the free energy. 

We now consider the case c = 0, c > 0. For this case, there is a one-parameter 
(a) family of lowest-energy states. Indeed, consider the following one-parameter family of 
densities: 

pi = as , pi = (1 — a)s , p2 = s + (1 — a)s * s , < a < 1 , (30) 

and all other particle and hole densities are equal to zero. These densities obey the 
constraints (7), and give (independently of the value of a) the same lowest value for 
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the energy. Moreover, to leading order in N, the spin is S z = 0, and the entropy is 
S = — ^ [a In a + (1 — a) ln(l — a)]. In particular, for a 7^ 0,1, the entropy is nonzero 
and is proportional to N, implying an infinite degeneracy of states. This degeneracy is 
consistent with the fact that, above the a = 1 vacuum, there are excitations (namely, holes 
in the sea of real roots) which have zero energy and non-zero momentum 2 . We are not 
aware of any other model with such properties 12 . 

We speculate that the system can be brought to these various states by first preparing 
the system at finite T and H, and then approaching the origin (T = , H = 0) of the 
(T , i7)-plane from appropriate directions. The ground state of the system is the state 
which is reached by approaching the origin along the line H = 0. Unfortunately, we 
cannot determine the particular value of a corresponding to this state, since this would 
entail computing liniT^o lirnff^o ei/T, while we know how to calculate only for T << H 
13 . By approaching the origin of the (T,i7)-plane along the line T = 0, the state with 
a = is reached. 

We have calculated the free energy for small values of T and H, with T « H. The 
calculation is similar to the one above, except that now one must take into account that 
£i(A) does not have a zero. We find that, to leading order, the free energy per site is given 
by 

f = e + -^H 2 \nH+ -^T 2 InH . (31) 
Atk a c 12tvc 

Contrary to appearance, this result does not imply that / diverges for H — > 0, since in 
the region where the calculation is valid (T << H), the last term is finite. This result is 
nevertheless unusual, since it implies that Ch/T is a function of H. 

Similar results are obtained for the case c = 0, c > 0. For the cases c < 0, c > and 
c > 0, c < , the model either has a ferromagnetic ground state and a finite gap, or it has 
no gap, depending on the precise values of c , c, and H. For the cases c, c < , the model 



is ferromagnetic. 



We thank L. Susskind for a very helpful discussion and F. Zuo for bringing Ref. 1 to 
our attention. One of us (H J de V) thanks the Department of Physics of the University 
of Miami for the warm hospitality extended to him. This work was supported in part by 
the National Science Foundation under Grant No. PHY-92 09978. 

Note added: After this work was completed, we received a preprint (Ref. 14) which 
discusses chains of alternating spin 1/2 and spin s. However, that paper considers neither 
the effect of an external magnetic field, nor the behavior of the models away from the 
conformally invariant point (i.e., for c^c). 
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